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SUMMARY  &  CONCLUSIONS 

Knowledge  aided  signal  processing  calls  for  the  use  of 
existing  information  such  as  digitized  image  maps  to  further 
suppress  clutter  in  a  typical  radar  detection  operation.  In  the 
context  of  space  based  radar  (SBR)  applications,  given  the 
SBR  parameters  and  the  look  angles  the  first  problem  is  to 
translate  these  parameters  into  a  data  point  on  a  map.  This  is 
realized  by  converting  these  parameters  into  conventional 
latitude  and  longitude  and  then  using  the  corresponding  image 
intensity  to  estimate  the  clutter  covariance  matrix.  Using  this 
approach,  it  is  possible  to  pre-compute  the  desired  weight 
vector  for  a  location  of  interest  on  earth  that  the  SBR  is 
assigned  to  visit.  As  a  result,  using  the  measured  data  together 
with  the  weight  vector,  it  is  possible  to  test  for  target  detection 
hypothesis  accurately  at  the  time  of  visit  by  the  SBR. 

1.  INTRODUCTION 

An  interesting  problem  in  knowledge  aided  Space  Time 
Adaptive  Processing  (KA-STAP)  is  how  to  make  use  of 
traditional  image  maps  to  cancel  out  clutter  and  interference 
components  while  looking  for  a  target.  This  approach  has 
applications  in  any  ‘search  and  detect’  operation  including 
airborne  and  space  based  radar  (SBR).  In  what  follows  we 
shall  illustrate  our  approach  for  the  SBR  applications.  The 
problem  in  this  context  is  to  translate  the  local  geometry  to  the 
SBR  platform  so  that  the  map  intensities  can  be  translated  into 
meaningful  clutter  covariance  estimates. 

In  its  general  setting  two  points  with  respective  latitudes 
and  longitudes  and  (a2,j82)  are  chosen  on  a  map,  the 

first  indicating  the  location  of  the  transmitter/receiver  such  as 
an  SBR  and  the  second  indicating  the  point  of  interest  (see 
Fig.  1).  In  addition,  given  the  SBR  height  H  and  its  orbit 

inclination  angle  77.  at  the  equator,  the  first  step  is  to 
determine  all  other  desired  parameters  such  as  SBR  elevation 
angle  0EL ,  SBR  array  azimuth  angle  6 Az  ,  the  true  ground 
range  R  ,  grazing  angle  y/ ,  the  array  footprint  on  earth,  range- 
foldover  points  RrR2,...  on  the  ground,  and  finally  the 

correction  factor  on  the  azimuth  due  to  the  earth’s  rotation 
(crab  angle)  (see  Fig.  2  and  Fig.  3)  [1,  2,  3]. 
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Fig.  1:  SBR  and  point  of  interest  co-ordinates 
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Fig.  2:  SBR-Earth  Geometry 
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Fig.  3:  SBR  2-D  Array  Gain  Pattern  and  Range  Foldovers 

The  SBR  Parameters  are  listed  below  (see  [4]  for 
derivations): 

R  =  Re  cos  1  (sin  ax  sin  a2  +  cos  a,  cos  a2  cos(/?2  — /?, ))  ( 1 ) 
oc  =  9e  =  R/Re  (2) 

Rs  =  A2  +  (Re  +  #  )2  -  +  # )  cos(fl  /  ^ )  (3) 


formulas  when  the  parameters  are  appropriately  incremented. 
Such  is  the  case  when  computing  the  clutter  covariance  matrix 
Rc  using  an  image  map. 

In  general  the  clutter  covariance  matrix  has  the  form  [1,2, 

3], 

Rc  =ZPkG(ik)s(ik,codt)SH(ik,a)dt)  (9) 

k 

where  Pk  is  the  intensity  of  the  kth  clutter  patch,  G(6k)  the 
corresponding  array  gain,  and  s(Ok,G)d  )  is  the  space-time 
steering  vector  associated  with  the  kth  patch.  Here  0k  refers  to 
both  azimuth  angle  0AZk  and  the  elevation  angle  dEL  of  the 

point  of  interest,  and  in  practice  are  incremented  to  cover  all 
guard  cells  used  to  estimate  Rc  in  (9).  In  what  follows  inverse 
formulas  for  latitude  and  longitude,  given  0EL  and  0AZ  ,  are 
developed. 

2.  AZIMUTH  INCREMENT 

From  Fig.  4,  we  have 

0'az=0az+A0az,  (10) 

where  0AZ  and  0AZ  refer  to  the  azimuth  angles  corresponding 
to  the  point  of  interest  before  and  after  the  increment, 
respectively. 


sin(^'/2  +  ^r)  _  sin<9e  _  sin  <9^ 
(Re+H)  ~~R^~  Re 


0EL  =  sin" 


K 

yRs 


sin  (R/  R) 


(4) 

(5) 


sin  y=  cos«2 


sin(A-A) 

sin  a 


d^  =  r+s-7z, 


(6) 

(7) 


Fig.  4:  Azimuth  increment 

Note  that  in  this  case  the  range  remains  the  same  and 
there  is  a  change  in  the  azimuth  angle  only.  Thus,  (7)  can  be 
used  to  compute  the  new  y . 


where 


=  sin 


cos  rji 


(8) 


These  formulas  give  the  SBR  parameters  in  terms  of 
(or, , /( ) ,  (or2,/?2),  and  the  SBR  height  H.  Using  these 
parameters  the  SBR  array  pattern  on  the  ground  can  be  plotted 
as  shown  in  Fig.  3.  Clearly  when  pointed  at  (or2,/?2)  the  SBR 
receives  data  from  all  points  that  are  illuminated  in  Fig.  3,  and 
hence  it  is  necessary  to  make  angular  increments  both  in  the 
elevation  and  azimuth  directions  to  access  the  neighbouring 
points.  As  a  result,  it  is  necessary  to  develop  the  inverse 


Fig.  5:  SBR-point  of  interest  spherical  geometry 
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In  the  spherical  triangle  ABC  (Fig.  5),  two  sides  and  an 
included  angle  y  are  known.  This  gives  [5] 

cosc  =  sin  ax  cos  a  +  cos  ax  sin  or  cos  y  (11) 

and  hence  the  desired  latitude  of  the  new  point  is  given  by 

at  =7i  12- c. 


the  mainbeam  at  the  point  of  interest,  as  shown  in  Fig.  7.  The 
associated  range  points  corresponding  to  these  returns 
constitute  the  range  ambiguities  for  the  point  of  interest. 

4.  RANGE  FOLDOVER 

The  ground  plane  resolution  of  the  SBR  is  limited  by  the 
pulse  frequency  and  the  grazing  angle  y/.  This  causes  range 
ambiguities  spaced  at  the  limiting  resolution  SR  ,  given  by 


Once  at  is  known,  /?.  can  be  obtained  from  the  relation  [5] 


~f\)  _  Sln  y 

sin  a  sin  c 


(13) 


3.  ELEVATION  INCREMENT 

In  this  case  the  new  elevation  angle  corresponding  to  the 
point  of  interest  is  given  by  (see  Fig.  6) 


O'tL  =®EL  +A6EL’ 


(14) 


where  0EL  refers  to  the  elevation  angle  before  the  increment. 


Fig.  6:  Elevation  increment 

With  the  incremented  value  of  the  elevation  angle,  Eq.  (4) 
can  be  used  to  determine  the  grazing  angle  yr  .  Thus, 


SR=—secy/.  (18) 

From  (18)  and  for  short  range  regions  where  the  grazing 
angle  y/  is  closer  to  ;r/2,  the  range  resolution  is  very  poor 
since  SR  —>  °° .  For  long  ranges,  the  resolution  approaches  its 
limiting  value  cT  12  as  y/  — >  0. 
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Fig.  7:  Range  Foldover 


y/  =  cos 


sin  eEL 

ll+tf/tfj 


and  the  relation 


If  Rs  represents  the  slant  range  at  the  end  of  one  pulse, 
(15)  then  Rs  +cT  /  2  is  the  new  slant  range  at  the  end  of  the  next 
pulse.  Let  R  +  AR  represent  the  new  range  corresponding  to 
the  second  pulse  shadow  on  the  ground.  Using  Fig.  2,  we 
have 


0e  =  nl2-6n  -y/  (16) 

can  be  used  to  determine  the  new  range  R  .  This  gives 

R  =  Re0e=Re(7T/2-0EL-ys).  (17) 

Note  that  0AZ  and  y  remain  unchanged  and  (11)  -  (13) 
can  be  used  to  compute  (tz.,/?.) .  Knowing  the  latitude  and 

longitude,  the  power  level  can  be  extracted  from  the  map.  This 
procedure  is  repeated  to  cover  all  guard  cells  of  interest  in  (9). 

Another  refinement  in  modelling  the  clutter  covariance 
matrix  is  to  introduce  the  effect  of  range  foldover  ambiguities 
into  (9).  Range  foldover  occurs  when  clutter  returns  from 
previously  transmitted  pulses  are  coherently  combined  with 


or 


(R+cTI  2f  = 

R2  +(Re  +H)2  -  2 Re (Re  +  H) cos 


fR+±y 

V  Re  J 


(19) 


=  Re  cos  1 


R]  +  (Re  +  Hf  -  (Rs  +  cT  /  2 )2 
""  2  R(R+H) 


-R.  (20) 


Note  that  A^  is  a  decreasing  function  of  R.  Moreover 
when  R  is  relatively  small,  the  distance  between  the  pulse 
shadows  on  the  ground  is  large  and  it  decreases  as  R  increases. 
However  for  large  values  of  range,  AR  remains  constant  at  its 
limiting  value  cTr! 2  .  As  a  result  and  as  the  range  increases, 
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this  has  the  additional  effect  of  packing  several  range 
ambiguities  within  a  mainbeam  footprint. 

If  we  assume  to  be  constant  within  a  mainbeam 
footprint,  then  the  number  of  range  ambiguities  within  the 
mainbeam  is  given  by 

Na=(RT-RH)/AR.  (21) 

Otherwise,  starting  at  the  heel  of  the  mainbeam  located  at 
range  RH  ,  ,  •  •  • ,  are  sequentially  calculated  with  the 

mainbeam  footprint  until  the  toe  is  reached.  Thus,  the  number 
of  range  ambiguities  Na  within  a  mainbeam  satisfies  the 
equation 

A rx  +  Ar^  h  I-  ARffa  —  RT  —  Rh  .  (22) 

Each  range  ambiguity  contribution  to  the  clutter 
covariance  matrix  has  the  same  form  as  in  (9).  As  a  result,  the 
clutter  covariance  matrix  that  accounts  for  range  ambiguities 
has  the  following  form 

K  =  I2X  G{9km)s(dkm,codim  )sH  {O^co^ ) .  (23) 

m- 0  k 

Here  Pkm  corresponds  to  the  intensity  of  the  kth  patch 

from  the  mth  range  foldover.  For  each  range  ambiguity  in  (23), 
the  azimuth  increment  method  described  earlier  can  be  applied 
to  compute  its  clutter  covariance  matrix  contribution. 

Finally  the  clutter  Doppler  o)d  in  (23)  can  be  modelled 

with  and  without  taking  earth’s  rotation  into  account.  The 
effect  of  earth’s  rotation  tends  to  modify  the  clutter  Doppler  in 
terms  of  a  crab  angle  (/)c  and  a  crab  magnitude  pc  affecting  its 

phase  and  magnitude  [2,  3].  With  T  representing  the  transmit 
pulse  period  and  Vp  the  satellite  velocity,  it  can  be  shown  that 
the  clutter  Doppler  with  earth’s  rotation  for  the  kth  patch  is 
given  by  [2,  3] 

2VT 

=  -  JJ -  Pc  sm  eEU  cos  (0AZm+</>c).  (24) 

5.  TERRAIN  MAPPING 

Using  the  method  described  above  to  update  the  latitude 
and  longitude,  terrain  data  from  any  stand  alone  image  map 
can  be  used  to  custom  compute  the  covariance  matrix.  For 
illustration  purposes,  we  have  used  image  intensity  data  from 
a  world  map  at  http://www.cs.brown.edu/exploratories. 

Fig.  8  shows  the  SBR  coverage  for  two  regions,  (i)  the 
eastern  United  States  (6^=34°,  ^=-100°)  and  (ii)  the 
Middle  East.  To  compute  the  clutter  covariance  matrix  for 
each  region,  the  array  patterns  and  the  range  foldovers  are 
projected  on  the  surface  of  the  earth  using  (10)-(17).  In 
addition,  the  range  foldovers  are  included  in  this  case  since  we 
are  using  a  pulse  repetition  frequency  of  500  Hz.  For  an  SBR 


at  a  height  of  506  km,  we  have  Rmax~2400km  and  Na  is 
equal  to  6  or  7  depending  on  the  target  range.  By  carrying  out 
the  latitude/longitude  -  azimuth/elevation  conversion,  the 
terrain  of  interest  can  be  superimposed  on  an  image  map  to 
obtain  the  true  terrain  image  intensities.  These  intensity  values 
Pik  together  with  the  terrain  geometry  described  earlier  can  be 
used  to  estimate  the  clutter  covariance  matrix  Rc  for  the  point 
of  interest  with  elevation  and  azimuth  angles  0ELq  and  GAZq  . 
Knowing  the  SBR  array/pulse  configuration,  if  s(Go,cod) 

represents  the  space-time  steering  vector  for  the  target 
location,  then 

w=R-1s(eo9G)d )  (25) 

represents  the  desired  weight  vector  that  can  be  applied  to  a 
data  sample  to  determine  the  presence  or  absence  of  a  target  of 
interest. 

Fig.  9  shows  the  performance  of  the  optimum  weight 
vector  in  terms  of  its  ability  to  perform  clutter  suppression 
with  and  without  the  presence  of  range  foldover  using  the 
Signal  to  Interference  plus  Noise  (SINR)  measure,  given  by 

SINR  =  sH (Go,0)d) R -1  s(0o,ct)d) .  (26) 

Their  combined  effect  results  in  appreciable  performance 
degradation  as  shown  in  Fig.  9.  In  this  context,  the  presented 
approach  gives  a  tool  to  evaluate  the  performance  degradation 
at  various  points  on  earth.  A  preliminary  simulation  shows  the 
clutter  suppression  for  different  points  of  interest  across  the 
United  States. 


Fig.  8:  SBR  coverage  and  range  foldover  maps.  Attached 
images  are  used  to  estimate  the  clutter  covariance  matrix. 
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SI  NR  ys  Velocity 


Fig.  9:  MDV  performance  with  and  without  Range  Foldover 
and  Crab  Angle  (Range  =  1200km,  Azimuth  =  88°) 


SI  NR  vs  Velocity 


Fig.  10:  Minimum  Detectable  Velocity  for  -5  dB  threshold. 


6.  EFFECTS  OF  TERRAIN  ON  CLUTTER  SIMULATION 

As  previously  described,  the  intensities  obtained  from  an 
image  map  can  be  used  when  estimating  the  clutter  covariance 
matrix,  R  .  To  test  these  effects,  multiple  simulations  were 
run  at  different  points  over  the  United  States  within  the  three 
regions  of  interest  as  shown  in  Fig.  11.  First,  the  target  range 
and  azimuth  angle  were  held  constant  at  900  km  and  88°  and 
the  location  of  the  SBR  was  moved  among  the  three 
longitudes  (-125°,  -110°,  -95°)  and  four  latitudes  (36°,  38°, 
40°,  42°).  By  making  use  of  (23)-(26)  the  clutter  suppression 
performance  at  each  of  these  locations  can  be  studied  using 
the  output  SINR.  The  results  are  shown  in  terms  of  the 
Minimum  Detectable  Velocity  (MDV)  measure  for  ground 


moving  target  detection.  From  Fig.  10  MDV  is  a  measure  of 
the  ability  to  detect  slow  moving  targets  in  the  presence  of 
clutter  and  noise.  From  there,  clearly  a  sharp  clutter  ridge 
provides  superior  performance.  For  the  sake  of  comparison, 
we  have  set  the  threshold  for  MDV  at  5  dB  below  the  noise 
floor  as  shown  in  Figs.  9-10.  Notice  that  when  earth’s  rotation 
is  present  MDV  is  undefined  as  shown  in  Fig.  9,  indicating  the 
a  severe  degradation  in  performance. 

The  average  MDV  values  for  the  three  regions  are  shown 
in  Fig.  11.  The  MDV  variations  at  different  latitudes  within 
each  region  are  shown  in  Fig.  12.  Notice  that  the  Midwest 
region  gave  an  average  MDV  of  23.2749  m/s  whereas  the 
West  region  gave  an  average  MDV  of  24.3179  m/s  (See  Fig. 
11).  The  East  region  MDV  was  also  slightly  higher  than  the 
Midwest.  The  Midwest  is  generally  flatter  and  more  uniform 
than  the  other  regions  so  one  would  expect  lower  MDV’s  due 
to  the  “homogeneous”  nature  of  the  background  terrain. 


Region  1 

Region  2 

Region  3 

(West) 

(Midwest) 

(East) 

MDV  = 

MDV  = 

MDV  = 

24.3179  mis 

23.2749  mis 

23.6780  mis 

SBR 

SBR 

SBR 

longitude  = 

longitude  = 

longitude  = 

-125° 

-110° 

-95° 

Fig  11:  Regions  used  for  Terrain  Mapping. 

Similarly,  we  repeated  the  above  simulations  while 
keeping  the  latitude  constant  at  40°  and  adjusting  the  azimuth 
around  90°  (See  Fig.  13).  In  this  case,  MDV  increase  above  the 
other  regions  for  the  Midwest  at  larger  azimuth  angles.  This 
could  be  due  to  terrain  variations  such  as  the  southern  end  of 
the  Rocky  Mountains  and  the  Gulf  of  Mexico.  Yet  as  before, 
MDV  values  for  the  Midwest  are  lower  when  the  beam 
covered  the  plains. 


Fig.  12:  Range  =  800km,  Azimuth  =  88° 
MDV  calculated  at  four  latitudes  in  each  region. 


0-7803-8882-8/05/$20.00  (C)  2005  IEEE 


REFERENCES 


Azimuth 


Fig.  13:  Range  =  800km,  latitude  =  40°.  MDV 
calculated  for  different  azimuth  in  each  region. 


7.  CONCLUSIONS 

This  paper  discusses  methods  for  incorporating  existing 
digital  maps  using  SBR-earth  geometry  to  generate  realistic 
clutter  covariance  nulling  schemes  for  ground  moving  target 
detection  applications.  SBR-earth  geometry  parameters  need 
to  be  updated  into  the  latitude-longitude  scale  to  determine  the 
actual  image  intensities.  This  approach  allows  pre-computed 
optimum  weight  vectors  to  be  determined  that  can  be  used  on 
a  single  snapshot  of  measurement  to  determine  the  presence  or 
absence  of  a  target.  Performance  results  using  this  approach 
across  the  United  States  terrain  are  reported  in  terms  of 
Minimum  Detectable  Velocity  (MDV)  and  indicate  the  clutter 
suppression  efficiency. 
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